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, Abstract 

> . 

d , An L-basis associated to a linear second-order ordinary differential operator L is an infinite 

sequence of functions {(^fc}^ sucn that Lift = for k = 0, 1, Lipk = k(k — l)ipt-2, for k — 2, 3, . . . 
and all ipk satisfy certain prescribed initial conditions. We study the transmutation operators 
related to L in terms of the transformation of powers of the independent variable {(a; — xo) k } k _ 
to the elements of the L-basis and establish a precise form of the transmutation operator realizing 
this transformation. We use this transmutation operator to establish a completeness of an infinite 
system of solutions of the stationary Schrodinger equation from a certain class. The system of 
solutions is obtained as an application of the theory of bicomplex pseudoanalytic functions and 
its completeness was a long sought result. Its use for constructing reproducing kernels and solving 
boundary and eigenvalue problems has been considered even without the required completeness 
ON ■ justification. The obtained result on the completeness opens the way for further development 

and application of the tools of pseudoanalytic function theory. 



1 Introduction 



Transmutation operators are a widely used tool in the theory of linear differential equations (see, e.g., 
PQ, [5], [22], [23], [35J and the recent review [H]). It is well known that under certain quite general 
conditions the transmutation operator transmuting the operator A = + q{x) into B = —j^z 

is a Volterra integral operator with good properties. Its kernel can be obtained as a solution of the 
Goursat problem for the Klein-Gordon equation with the variable coefficient. In the book [IT] another 
approach to the transmutation was developed. It was shown that to every (regular) linear second- 
order ordinary differential operator L one can associate a linear space spanned on a so-called L-basis 
- an infinite family of functions {^fel^Lo sucn that Lifk = for k = 0, 1, Lipk = k(k — l)(fk-2, for 
k = 2, 3, . . . and all tpk satisfy certain prescribed initial conditions. Then the operator of transmutation 
was introduced as an operation transforming functions from one such linear space corresponding to 
a certain operator L to functions from another linear space corresponding to another operator M, 
and the transformation consists in substituting the L-basis with the M-basis preserving the same 
coefficients in the expansion. 

In the present work we find out how the canonical Volterra integral transmutation operator acts 

, ,2 

on powers of the independent variable x (which represent a basis associated with the operator -t-tt), 
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introduce a parametrized family of transmutation operators and construct a transmutation operator 
which transforms the powers x k into the functions fk from the L-basis. We prove that it is indeed a 
transmutation and can be written in the form of a Volterra integral operator. 

We apply this result to prove the completeness of certain families of solutions of linear two- 
dimensional elliptic equations with variable complex-valued coefficients. These families of solutions 
were obtained earlier [7], [T7] as scalar parts of bicomplex pseudoanalytic formal powers and used 
for solving boundary value and eigenvalue problems [5]. Nevertheless no result on their complete- 
ness even in simplest cases was known due to profound differences between complex and bicomplex 
pseudoanalytic function theories and inapplicability of many classical results and techniques in the 
bicomplex situation. The use of the constructed transmutation operators and their extremely fortu- 
nate transformation properties regarding the L-bases allow us to observe that the infinite families of 
solutions mentioned above are nothing but the transmuted harmonic polynomials. Using their well 
known completeness properties together with the properties of the Volterra integral transmutation 
operators we obtained several results on the completeness of families of solutions for equations with 
variable complex-valued coefficients. 

In Section 2 we introduce the definition and some basic facts about transmutations together with 
an example which we constructed for illustrating some results of the present work. In Section 3 we 
introduce the L-basis as a system of recursive integrals. In Section 4 we study the action of the 
transmutation operators on the recursive integrals and construct the transmutation operator which 
transforms powers of x into functions of the L-basis. In Section 5 we introduce several definitions 
and results from the recently developed bicomplex pseudoanalytic function theory and explain its 
relation to linear second-order elliptic equations with variable complex-valued coefficients. In Section 
6 we construct infinite families of solutions for a class of such equations and show that they are 
images of harmonic polynomials under the transmutation operator. We use this fact to prove their 
completeness under certain additional conditions. Section 7 contains some concluding remarks. 

2 Transmutation operators for Sturm-Liouville equations 

According to the definition given by Levitan (22], let E be a linear topological space, A and B be 
linear operators: E — > E. Let E\ and E2 be closed subspaces of E. 

Definition 1 A linear invertible operator T defined on the whole E and acting from E\ to E% is called 
a transmutation operator for the pair of operators A and B if it fulfills the following two conditions. 

1. Both the operator T and its inverse T^ 1 are continuous in E; 

2. The following operator equality is valid 



Our main interest concerns the situation when A — —-§^1 + q(x), B = —-^7, and q is a continuous 
complex- valued function. Hence for our purposes it will be sufficient to consider the functional space 
E = C 2 [a,b] with the topology of uniform convergency. For simplicity we will assume that the interval 
is symmetric with respect to the origin, thus E = C 2 [— a, a]. 

An operator of transmutation for such A and B can be realized in the form (see, e.g., [22] and 



AT = TB 



(1) 



or which is the same 



A = TBT 



1 



23) 




(2) 



where K{x, t) is a unique solution of the Goursat problem 
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/ d 2 \ d 2 

— n(r\ K t. t\ = — 



y dx 2 q(x))K(x,t) = w K{x,t), (3) 

1 r 

K(x,x) = - q(s)ds, K(x,-x)=0. (4) 

2 Jo 

An important property of this transmutation operator consists in the way how it maps solutions 
of the equation 

v" + uj 2 v = (5) 

into solutions of the equation 

u" — q(x)u + lj 2 u = (6) 
where w is a complex number. Denote by eo(icu,x) the solution of ([6]) satisfying the initial conditions 

eo{iu>, 0) = 1 and e (iu>,0) = iu. (7) 

The subindex "0" indicates that the initial conditions correspond to the point x = and the letter 
"e" reminds us that the initial values coincide with the initial values of the function e lux . 
The transmutation operator ([2]) maps e lu)x into eo(iuj,x), 

e (iuj,x) = T[e iux ] (8) 

(see [23, Theorem 1.2.1]). 

Following [23] we introduce the notations 

K(x, t;h) = h + K(x, t) + K(x, -t)+h J {K(x, £) - K(x, ~0} d t 

where h is a complex number, and 

K(x, t; oo) = K(x, t) - K(x, -t). 

Theorem 2 f23j/ Solutions c(w, x; h) and s{u, x; oo) of equation satisfying the initial conditions 

c(w, 0; h) = 1, c' x (w, 0;h) = h 

s(w, 0; oo) = 0, s^(w,0;oo) = l 

can be represented in the form 

c[u), x; h) = cosojx + / K(x, t] h) coswi dt (9) 
Jo 

and 



, . sinwx f . sincjt 
s(w,x;oo) = V \ K{x,t;oo) dt. (10) 

UJ In LO 



The operators 



and 



T c u(x) — u(x) + / K(x,t;h)u(t)dt 
Jo 

T s u(x) = u(x) + / K(x,t;oo)u(t)dt 
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are not transmutations on the whole space C 2 [~a,a], they even do not map all solutions of (J5J) into 
solutions of ([6]). For example, as we show below 



■2 



+'/(•'•) )7M] = 7\ 



d 2 
■ — (I) 

dx 2 1 ' 



= 



when q is constant. 



,2 j2 

Example 3 Transmutation operator for operators A := -P-^ + c, c is a constant, and B := -rr. 
According to ]23\. (1.2.25), (1.2.26)], finding the kernel of transmutation operator is equivalent to 
finding the function H(s, t) = K(s +t,s — t), satisfying Goursat problem 

d 2 H{s,t) , , cs , 

d K sdt ' =-cH(s,t), H(8,0) = --, H(0,t) = 0. 

The solution of this problem is given by \ 121 (4-85)] 

H(s,t) = -\ [ J^^T)) - - ^ J f^\ 

where Jq and J± are Bessel functions of first kind, and the formula is valid even if the radicand is 
negative. Hence, 



1 2 ' 2 / 2 x-y ' 



From (jllj) we get the 'sine' kernel 



tJc{x 2 -t^MJcix 2 ~t 2 )) , s 

K (x,t;oo 12 

x z — t A 



and can check the above statement about operator T s , 



T a [lJ(x) = 1- / 2~72 -dt= J {xy/c), 



t^c{x 2 ~t 2 ).h{^c{x 2 -t 2 )) 
x 2 -t 

, „Uhl_ y/cJi{xy/c) 

dx 



3 A complete system of recursive integrals 

Let / e C 2 (a, b) n C 1 [a, b] be a complex valued function and f(x) ^ for any x g [a, b]. The interval 
(a, 6) is supposed to be finite. Let us consider the following auxiliary functions 

X^(x) = X^(x) = 1, (13) 

a; 

xW(x) =n[ X^Hs) (f 2 (s))^ ir " da, (14) 
X^(x) = nlx^- 1 \s)(f(s)) { - ir da, (15) 



where x$ is an arbitrary fixed point in [a, b]. We introduce the infinite system of functions {</?fc}fcLo 
defined as follows 
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f(x)X^(x) 7 fcodd, 
f(x)X^(x), k even, 



where the dchnition of X^> and X^ is given by (p~3]l - (p~5|) with xo being an arbitrary point of the 
interval [a, 6]. 

Example 4 Let f = 1, a = 0, b = 1. TTien it is easy to see t/iat choosing xo = we have ifk{x) = x k . 
k € No where by No we denote the set of non-negative integers. 

In [TS] it was shown that the system {<Pfe}^L is complete in £2(0, &) and in [TO] its completeness 
in the space of piecewise differentiable functions with respect to the maximum norm was obtained 
and the corresponding series expansions in terms of the functions ifk were studied. 

The system (|T5)) is closely related to the notion of the L-basis introduced and studied in [TT| . Here 
the letter L corresponds to a linear ordinary differential operator. This becomes more transparent 
from the following result obtained in [16] (for additional details and simpler proof see [17] and [20] ) 
establishing the relation of the system of functions {vfcj^o ^° Sturm-Liouville equations. 

Theorem 5 Let q be a continuous complex valued function of an independent real variable 

x G [a, 6], A be an arbitrary complex number. Suppose there exists a solution f of the equation 

f" -qf = (17) 
on (0,6) such that f £ C 2 [a, b] and f 7^ on [a,b]. Then the general solution of the equation 

U — qu = Xu (18) 

on (a, b) has the form 

u = ciui + c 2 u 2 
where C\ and c 2 are arbitrary complex constants, 

00 yk 00 yk 

Ul = E(2fc)!^ and U2 =£^^<^ +1 (19) 
and both series converge uniformly on [a, 6]. 

Remark 6 Lt is easy to see that by definition the solutions U\ and u 2 satisfy the following initial 
conditions 

ui(x ) = f(x ), u[(xo) = f'(x ), (20) 
u 2 {x )=Q, u' 2 (x )=l/f(x ). (21) 



4 Transmutations and systems of recursive integrals 

Let us obtain the expansion of the solution eo(iw, x) from Section[5]in terms of the functions ip-^. We 
suppose that / is a solution of (fTT)) fulfilling the condition of Theorem [5] on a finite interval (—a, a). 
We normalize / in such a way that /(0) = 1 and let /'(0) = h where h is some complex constant. 
Then according to Remark [5] the solutions (fT9"|) of equation (TT8"|) have the following initial values 

m(0) = 1, «i(0) = h, tta(0) - 0, u' 2 (0) = 1. 

Hence due to ([7]) we obtain eo(iu>,x) = u 1 (x) + (iuj — h)u 2 (x). From ([SJ and (fT9| we have the equality 

k=0 y ' k=0 y ' j=0 J J-x \ -_ Q J J 
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As the series under the sign of integral converges uniformly and the kernel K(x,t) is at least contin- 
uously differentiable (for a continuous q [23] ) we obtain the following relation 



The equality holds for any lj hence we obtain the termwise relations 

ip k = T[x k ] when k is odd (22) 

and 

ifk — -, T^k+i = T\x k ] when k £ No is even. (23) 

k + 1 

Taking into account the first of these relations the second can be written also as follows 



fk = T 



k + l 



when fc e No is even. (24) 



Thus, we proved the following statement. 



Theorem 7 Let q be a continuous complex valued function of an independent real variable x £ 
[—a, a], and f be a particular solution of JJ7[ j such that f £ C 2 [— a, a], / ^ on [—a, a] and normalized 
as /(0) = 1. Denote h := /'(0) G C. Suppose T is the operator defined by |IJ) where the kernel K 
is a solution of the problem |3p, and ifk, k £ No are functions defined by U6\) . Then equalities 
hold. 



Thus, we clarified what is the result of application of the transmutation T to the powers of the 
independent variable. This is very useful due to the fact that as a rule the construction of the 
kernel K(x,t) in a more or less explicit form up to now is impossible. Our result gives an algorithm 
for transmuting functions which can be represented or at least approximated by finite or infinite 
polynomials in the situation when K(x,t) is unknown. 



Remark 8 Let f be the solution of satisfying the initial conditions 

/(0) = 1, and /'(0) = 0. (25) 

If it does not vanish on [—a, a] then from Theorem^ we obtain that fk — T\x k \ for any k € No- 
In general, of course there is no guaranty that the solution with such initial values have no zeros on 
[—a, a] and hence the operator T transmutes the powers of x into ifk whose construction is based on 
the solution f satisfying \25)) only in some neighborhood of the origin. 

Similarly to Theorem [7] we obtain the following statement. 

Theorem 9 Under the conditions of Theorem^ the following equalities are valid 

<fik = T c [x k ] when k G No is even (26) 

and 

ip k =T s [x k ] when k £ N is odd. (27) 

Proof. It is easy to see that c(u),x;h) = u\{x) and s(w, x;oo) = u%(x) where u\ and u-i are defined 
by (|I9p . From here and from ([9]), (jlOp by expanding coswx and sinwa; into their Taylor series we 
obtain flU and ([27]). ■ 

Now, for a given nonvanishing solution of (|17|) on (—a, a) satisfying the initial conditions f(0) = 1 
and /'(0) = h where h £ C and for the corresponding system of functions (|16p we construct a 



G 



transmutation operator for the pair 4^ and — q{x) such that x k are transformed into ipk{x) on 
the whole segment [—a, a] for any k £ No- For this we introduce the following projectors acting on 
any continuous function (defined on [—a, a]) according to the rules P e f(x) — (f(x) + /(— x))/2 and 
Pof{%) — (f( x ) ~ f{~ x ))/2- Consider the following operator 



It is easy to see that by construction for an even k we obtain T[x k ] = T c P e [x k ] — T c [x k ] — iff. and 
analogously for an odd k, T[x k ] — ipk due to (|27p . Thus, tpk = T[x k ) for any k £ No. Moreover, the 
operator T can be written as a Volterra operator in a form similar to @ . We have 



Tu(x)=u(x) + K(x,t;h)u(t)dt (28) 

where 

K(x, t;h) = ~ + K{x, t) + - [ (K(x, s) - K{x, -s)) ds. (29) 
2 2 

Let us notice that K(ir, t; 0) ~ K(x, t) and that the expression 

h r* 

K(x, t; h) - K(x, -t; h) = K(x, t) - K(x, -t) + - / (K(x, s) - K(x, -s)) ds = K(x, t) - K(x, -t) 

2 J-t 

does not depend on h. Thus, we obtain a way to compute K(ir, t; h) for any h by a given K(a;, t; hi) 
for some particular value hi. 

Theorem 10 The integral kernels K(x, t;h) and K.(x,t;hi) are related by the expression 

K(x,t;h)= h ~ hl +K(x i t;hi)+ h ~ hl [ (K(i, s; hi) - K(x, -s; hi)) ds. (30) 
2 2 Jt 

Let us prove that T is indeed a transmutation. 

Theorem 11 Under the conditions of Theorem^let us assume additionally that q £ C l [—a, a]. Then 
the operator (28\) with the kernel defined by i29\) transforms x k into fk(x) for any k £ Nq and 



■2 



^2 +'/'■'•) ) T 



2 



—(u) 
dx 2 



(31) 



for any u £ C 2 [—a,a\. 



Proof. Under the condition q £ C 1 [— a, a], the kernel K (x, t) in @ is twice continuously differcntiablc 
with respect to both x and t [531 Theorem 1.2.2]. Hence, the kernel K(x, t; h) is also twice continuously 
differentiable with respect to both x and t, and C 2 [— a, a] is invariant under the operator T, that is, 
the left-hand side of (|3Tj) is well defined for all u £ C 2 [— a, a]. 

Since (-^ + q(x)j <p k = k{k - l)(p k -2, k > 2, and + q(xj) ^ k = 0, k = 0, 1 (see [H]), 

the equality ([3T|) is valid for all powers x k and, by linearity, for all polynomials. Let u £ C 2 [—a,a\. 
Then, u" £ C[—a,a] and by the Weierstrass theorem there exists a sequence of polynomials Q n such 
that 

Q n — > u , n — > oo uniformly on [—a, a]. (32) 
Integrating (l3"2"j) twice we conclude that the sequence of polynomials defined by 

P n (x) = u(0) + u'(0)x+ [ [ Q n (s)dsdt 

Jo Jo 
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is such that 

P„ — > u, P n — > u and P n — > u" , n — > oo 

uniformly in [—a, a]. Since the kernel K.(x,t;h) is twice continuously diffcrcntiablc, it is easy to see 
that also 

T[P n ] -> T[u], (T[P n ])' -> (T[«])' and (T[P„])" -> (T[u])", n -j. oo 
uniformly in [—a, a]. Therefore, 



+ q (x)) T M = Bm (-*- +ff(a5 )]T[7',,] In,, T 



<te 2 



d 2 



It is possible to give another final part of the proof, revealing some properties of the adjoint 
operator T* . 

Proof. Consider u e C 2 [— a, a] and let {p n }neN be a sequence of polynomials such that p n — >• u, n — > 



Consider /„ = ^—^7 + <zOc)l T[p n ]. Since p n is a polynomial, we also have /, 

Let / = ^— -j^j + q(x)j T[u], f = T — -j£x{u) . From now on consider u, p n , /„, /, / as elements 

of the Hilbert space Li\— a, a]. Since T is the Volterra operator, it is continuous in L2[— a, a] space. 
Hence, Tp„ — > Tu, n — > 00. Consider any function ip € C 2 [— a,a], that is twice continuously 
differentiable and supported on some [a, ft] C (—a, a). Then we have 



(/»-/,V0= ^-^2 J (Tp n -Tn),^J = ^Tp B -Ttt, ->0, n ^ oc. 

(33) 

For the right-hand side of (|3ip . we need the adjoint operator T*. Since T is the Volterra 
operator, its adjoint is given by the expression |131 Chapter 3, Example 3.17] T*u(x) = u{x) + 
f_] x| K(t, X] h)u{t)dt + Jj^, K(i, h)u{t)dt, from which it is easy to see that T*V> G C 2 [— a, a] for any 
tfi € Cg[— a,a]. Hence, 



(/»-/» 



dx 2 



Pn-tt,-^a (T» ) ->0, n->oo. (34) 



It follows from (|3"3")l and ([M]) that (f — f,ip) — for any ^ G Cq[— a, a]. Since the set C 2 [— a, a] is 
dense in Li\— a, a], we have / = / as elements of -L2[ — a, a] as well as continuous functions. ■ 

Example 12 Consider the same operators A and B as in Example [3J Then the transmutation 
operator T is defined by the integration kernel (|29|) 



K(x,t;h) 



1 ^c(x 2 -y 2 )J 1 (^/c(x 2 -y 2 )) h 



h 



x-y 



2 JoW<x 2 - V 2 )) = K(x,t) + -J (Vc(x 2 - V 2 )). 



If we consider the function f(x) = e lKX , k 2 = c as a solution of Af = satisfying f(0) = 1, 
/'(0) = in, then the first four functions <pk are 

. , j„„ . . sm(nx) . . Kxe tKX — sin(Kx) , . 3(sin(«;x) — nx cos(kx)) 
<p fa) = e iKX , iptix) = — <p2(x) = . = tp a (x) = - — i— 

and /rom Theorem we obtain the integrals 



ipk(x) = x - -J I ^— ' h«y J (k\/(x 2 - y 2 )) I dy 



which validity can be checked numerically. 



8 



5 Bicomplex numbers and pseudoanalytic functions 



Together with the imaginary unit i we consider another imaginary unit j, such that 

j 2 = i 2 = — 1 and ij = j i. (35) 

We have then two copies of the algebra of complex numbers: Q := {a + ib, {a, b} C K} and Cj := 
{a + jb, {a, 6} C M}. The expressions of the form w = u+jv where {u, v} C C, are called bicomplex 
numbers. The conjugation with respect to j we denote as follows w — u — jv. The components u and 
v will be called the scalar and the vector part of w respectively. We will use the notation u = Sc w 
and v = Vecw. 

The set of all bicomplex numbers with a natural operation of addition and with the multiplication 
defined by the laws (|35[) represents a commutative ring with unit. We denote it by B. It contains 
zero divisors: the nonzero elements w such that ww — 0. Introducing the pair of idempotents 
P+ = §(l + ij) and P~ = {(P^f = P ± ) it is easy to see (e.g., OH p. 154]) that w = u + jv 

is a zero divisor if and only if w = 2P + u or w — 2P~u. For other algebraic properties of bicomplex 
numbers we refer to [21], [25] . 

We consider IB-valued functions of two real variables x and y. Denote d = + 3tL) and 



d = — ■ ^ n equation of the form 



2 V dx ~ J dy ' 

dx •* dy ' 

dw = aw + bw, (36) 



where w, a and b are B- valued functions is called a bicomplex Vekua equation. When all the involved 
functions have their values in Cj only, equation (|31)|) becomes the well known complex Vekua equation 
(see [T7], [35]). We will assume that w € C 1 (fi) where Vl C M. 2 is an open domain and a, b are Holder 
continuous in VL. _ 

When a = and b = & where : H Q possesses Holder continuous partial derivatives in Q 

and 4>(x, y) ^ 0, V(x, y) £ O we will say that the bicomplex Vekua equation 

dw — —y-w (37) 

<P 

is a Vekua equation of the main type or the main Vekua equation. 

For classical complex Vekua equations Bers introduced [3] the notions of a generating pair, gen- 
erating sequence, formal powers and Taylor series in formal powers. As was shown in [7], [17] the 
definition of these notions can be extended onto the bicomplex situation. Here we briefly recall the 
main definitions. 

Definition 13 A pair oj 'B- valued functions F and G possessing Holder continuous partial derivatives 
in O with respect to the real variables x and y is said to be a generating pair if it satisfies the inequality 

Vec(FG) ^ m Q. (38) 

Condition (|38[) implies that every bicomplex function w defined in a subdomain of f2 admits the 
unique representation w = <pF + %pG where the functions <f> and ip are scalar (C^-valued) . 

Remark 14 When F = 1 and G = j the corresponding bicomplex Vekua equation is 

dw = 0, (39) 

and its study in fact reduces to the complex analytic function theory. This is due to the fact that 
the functions P + w and P~w are necessarily antiholomorphic and holomorphic respectively. Indeed, 
application of P + and P~ to \39\) gives us 

d z P + w = and c\P~w = (40) 
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where d z — \{-§^~i-§^) anddz — ^(^:+*^)- Moreover, P + w = P + (u+jv) = P + (u—iv) and P w = 
P~(u + iv). Due to |^fl| ) the scalar functions w + := u — iv andw~ := u + iv are antiholomorphic and 
holomorphic respectively. We stress that w + is not necessarily a complex conjugate of w~ (u and v 
are <Ci-valued) . 

In general a reduction of the bicomplex Vekua equation (|36[) to a pair of decoupled complex 
Vekua equations is impossible. Application of P + and P~ to (|36| reduces it to the following system 
of equations 

d z w + = a + w + + b + w~ 

and 

c\w~ — a~w~ + b~w + 

for two complex functions w + and w~ with complex coefficients a^, b^. 
Assume that (F, G) is a generating pair in a domain f2. 

Definition 15 Let the M-valued function w be defined in a neighborhood of zo £ C Cj. In a 
complete analogy with the complex case we say that at z$ the function w possesses the (F, G)-derivative 
w(zq) if the (finite) limit 

s .. w{z) - XqF(z) - n a G(z) 
w(zo) = hm (41) 

z->ZO Z — Zo 

exists where Ao and fi$ are the unique scalar constants such that w(zq) = XqF(zo) + /iqG(zq). 

Similarly to the complex case (see, e.g., [TT1 Chapter 2]) it is easy to show that if w(zq) exists 
then at Zq, dw and dw exist and equations 

dw = a^ F ,G) w + 6(f,g)W (42) 

and 

w = dw — A(pQ)W — B^p G yw (43) 

hold, where <X(f,g)> ^>{f,g)i ^{f,g) an d Brp,G) are the characteristic coefficients of the pair (F,G) 
defined by the formulas 

FdG-GdF , FdG-GdF 
b 



FG-FG y ' ' FG-FG 

. _ FdG-GdF FdG-GdF 

A(F > G) - ~ FG-FG ' B(F > G) ~ FG-FG ' 

Note that FG-FG = -2j Vec(FG) ^ 0. 

If dw and dw exist and are continuous in some neighborhood of Zq, and if (|42|) holds at zq, then 

w(zo) exists, and (14*31 holds. Let us notice that F and G possess (F, G)-derivatives, F = G = and 
the following equalities are valid which determine the characteristic coefficients uniquely 

dF = a,(F,G)F + b(F,G)Fi 9G = a^ FG ^G + b( FjG )G, 

dF = A( FG) F + B( F t G)F, dG = A^ FyG )G + B( FG) G. 

If the (F, G)-derivative of a B- valued function w — <fiF + ipG (where the functions <f> and ip are scalar) 
exists, besides the form (J43J) it can also be written as follows w = dcj) F + dip G. 

Definition 16 Let {F,G) and (iq,Gi) - be two generating pairs in VL. (iq,Gi) is called successor 
of (F, G) and (F, G) is called predecessor of (iq, Gi) if 

a (Fi,Gi) = a (F,G) an d b(p lt Gi) = —B(F,G). 



10 



By analogy with the complex case we have the following statement. 

Theorem 17 Let w be a bicomplex (F, G)-pseudoanalytic function and let (Fi, G\) be a successor of 
(F,G). Then w is a bicomplex (Fi, G\)-pseudoanalytic function. 

Definition 18 Let (F,G) be a generating pair. Its adjoint generating pair (F,G)* — (F*,G*) is 
defined by the formulas 

2F G , 2G 



FG-FG FG-FG 
The (F, G)-integral is defined as follows 

J Wd {FtG) z=^F(z 1 )Sc J G*Wdz + G( Zl )Sc J F*Wdz^j 

where T is a rectifiable curve leading from z to z\. 

If W = (f)F + ipG is a bicomplex (F, G)-pscudoanalytic function where <j> and i(> are complex valued 
functions then 

Wd {F , G) z = W(z)-<j>(z )F(z)-ij(z )G(z), (44) 



/ 



and as F —G = 0, this integral is path-independent and represents the (F, G)-antiderivative of W. 

Definition 19 A sequence of generating pairs {(F m , G m )}, m — 0, ±1, ±2, . . . , is called a generating 
sequence if (F m+ i, G m+ \) is a successor of (F m ,G m ). If (Fq,Go) = (F,G), we say that (F,G) is 
embedded in {(F m ,G m )}. 

Let If be a bicomplex (F, G)-pseudoanalytic function. Using a generating sequence in which 
(F, G) is embedded we can define the higher derivatives of W by the recursion formula 

W®=W; W l™+i] = d (FM Wlm \ m= i, 2 ,.... 

dz 

Definition 20 The formal power Z$ '(a, Zo; z) with center at zq e fl, coefficient a and exponent 
is defined as the linear combination of the generators F m , G m with scalar constant coefficients X, ji 
chosen so that XF m (z ) + ^iG m (z ) = a. The formal powers with exponents n = 0, 1, 2, . . . are defined 
by the recursion formula 

z£ +1 Ha,z ;z) = (n + l) f Z^ +1 (a,z ; C)rf(F m , Gm )C- (45) 

J z 

This definition implies the following properties. 

1. Zm\a, z a ; z) is an (F m , G m )-pseudoanalytic function of z. 

2. If a' and a" are scalar constants, then 

Z&V + ja", z ; z) = a'zW(l, z ; z) + a"Z%\j, z ; z). 

3. The formal powers satisfy the differential relations 

d(F m ,G m ) Z m\a,Z ;z) r?(n-l), \ 

j- z = nZ ) n +i («. z o; Z). 
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4. The asymptotic formulas 

Z mH a ' z 0', z ) ~ a ( z ~ z o) n , z^zq 

hold. 

The case of the main bicomplex Vekua equation is of a special interest due to the following relation 
with the stationary Schrodinger equation. 

Theorem 21 jl^l Let W — W\ + jW 2 be a solution of the main bicomplex Vekua equation 

dW = in Cl (46) 

where W\ — ScW, W 2 = VecW and the Ci-valued function <j) is a nonvanishing solution of the 
equation 

- Au + q x (x,y)u = in CI (47) 



where q\ is a continuous Cj -valued function. Then W\ is a solution of y 7| ) in CI and W 2 is a solution 
of the associated Schrodinger equation 



Av + q 2 (x,y)v = in CL (48) 



where q 2 = g - qi- 



We need the following notation. Let w be a B-valued function defined on a simply connected 
domain CI with wi — Sew and w 2 = Vec w such that 

^-7^=0, V(«,y)6n, (49) 

and let T C be a rectifiable curve leading from (xo,yo) to (x,y). Then the integral 

Aw(x,y) := 2 (J w\dx + w 2 dy 

is path- independent, and all Ci-valued solutions ip of the equation dip = w in Cl have the form 
ip(x,y) — Aw(x,y) + c where c is an arbitrary C^-constant. In other words the operator A denotes 
the well known operation for reconstructing the potential function from its gradient. 

Theorem 22 'I4I Let W\ be a Ci-valued solution of the Schrodinger equation J^7| ) in a simply 
connected domain Cl. Then a Ci-valued solution W 2 of the associated Schrodinger equation J^ffp such 
that W\ + jW 2 is a solution of J^6'| ) in Cl can be constructed according to the formula 



W 2 = -A^j<t> z d\^ 

where c\ is an arbitrary Ci-constant. 

Vice versa, given a solution W 2 of |^<S| ) ; the corresponding solution W\ of J^7| ) such that W\+jW 2 
is a solution of |^6| ) has the form 

Wi = -<pA-(-^d(<f>W 2 )\+c 2 <p 
where c 2 is an arbitrary Ci-constant. 
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As was shown in [T2] (see also [T7]) a generating sequence can be obtained in a closed form, for 
example, in the case when <j) has a separable form (p — S(s)T(t) where s and t are conjugate harmonic 
functions and S, T are arbitrary twice continuously differentiable functions. In practical terms this 
means that whenever the Schrodinger equation (|47l) admits a particular nonvanishing solution having 
the form (f> — /(£) g(r)) where (£, 77) is one of the encountered in physics orthogonal coordinate systems 
in the plane a generating sequence corresponding to pp| can be obtained explicitly [17, Sect. 4.8] . The 
knowledge of a generating sequence allows one to construct the formal powers following Definition 1201 
This construction is a simple algorithm which can be quite easily and efficiently realized numerically 
[5], [5]. Moreover, in the case of a complex main Vekua equation which in the notations admitted in 
the present paper corresponds to the case of <p being a real- valued function (then the main bicomplex 
Vekua equation decouples into two main complex Vekua equations) the completeness of the system 
of formal powers was proved [5] in the sense that any pseudoanalytic in Q and Holder continuous 
on dfl function can be approximated uniformly and arbitrarily closely by a finite linear combination 
of the formal powers. The real parts of the complex pseudoanalytic formal powers represent then a 
complete system of solutions of one Schrodinger equation meanwhile the imaginary parts give us a 
complete system of solutions of the associated Schrodinger equation. 

In the bicomplex case the system of formal powers is constructed in the same way as in the 
complex situation and the system of functions 

{ScZ£°(l, *>;*), ScZ^\j,zo;z)}°° (50) 

l J n—Q 

is an infinite system of solutions of (|4T| . Nevertheless up to now no result on the completeness of the 
system of bicomplex formal powers or of the family of solutions (|50p has been proved. The reason 
is that such important basic facts which are in the core of pseudoanalytic function theory as the 
similarity principle are not valid for bicomplex Vekua equations. 

In the present work in order to establish such completeness for a certain class of Schrodinger 
equations we implement the transmutation operators. 



6 Complete families of solutions 

First we consider an important special case. Assume that (j) has the form <j){x, y) = f(x)g(y) where / 
and g are arbitrary Ci-valued twice continuously differentiable and nonvanishing functions defined on 
the segments [ai, 61] and [02, ^2] respectively. In this case there exists a periodic generating sequence 
with a period two corresponding to the main Vekua equation (|46p : 

(F,G) = (fg,4-), (F 1 ,G 1 )=(9-, J -£), (F 2 ,G 2 ) = (F,G) , (F a , G 3 ) = (F u G ± ), . . . , 



19/ \J 9 

and the corresponding formal powers admit the following elegant representation [3] . We consider the 
formal powers with the centre at the point (a:o,yo) G [<Zi,&i] x [02,^2]- We assume that f(xo) — 
5(2/0) = 1 and define the recursive integrals according to (|T3")) - (fT5"|) as well as the system of functions 
{^fel^Lo according to (fTo| . In a similar way we define a system of functions {^fc}^L corresponding 
to g, 

, , s f 9(y)Y^(y), fcodd, 
I 9(y)Y ( k> (y), k even, 



where 



y(°)(y) = y(°)(y) = l, (52) 
v 

? W (V)=n[ V ( "- 1} (s) (g 2 ( S )) ( - ir_1 ds, (53) 
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Y {n \y) =n /r ( "- 1) (s)( 5 2 (s)) ( 1} " ds. (54) 



Then the formal powers corresponding to (j46[) can be defined as follows. For a = a' + ia" and 
zq = xq + jyo we have 

2 (n) (a, *>,*) = /(x)ff(y) Sc .ZW ( a , ZQj z ) + (a;) J ( ^ Vec (a, z , z) (55) 

where 



t Z^(a,z ,z) = 

fe=o ^ ' 



(56) 



+ i a "Yl ( I ) x( "" fc) i kYik) for an odd 

fc=0 



and 



fc=0 



,Z (n) (a, 2„, z) = <T L ) jfey(fe) (57) 



ja" 



k=0 



^r W(n-^yW for an even 



Remark 23 Formulae \55\) - [57\ ) clearly generalize the binomial representation for the analytic pow- 
ers a(z — Zo) n . If one chooses f = 1 and g = 1 then Z^(a, Zq, z) — a(z — Zo) n . 



Consider the family of solutions of (147)) obtained from the scalar parts of the formal powers ([50 
We have 



Sc (1, z Q ; z) = <j>(x, y) Sc £ (f) X^j k Y^ 

k=0 ^ ' 



= f{x)g(y) J2 [ k \x {n - k) j k Y {k) for an odd i 

even k— 

ScZW(l,^;«) = /(a ! ) ff (y) ]T for an even 



A:=0 



ScZ< n \j,z ;z) = f(x)g(y) £ x("- fe )j fe+1 Y« for an odd 



odd k=l 



ScZ (n) (j,z ;z) =f(x)g(y) [ ) X {n - k) j k+1 Y {k) for an even n. 



odd fc=l v 

Taking into account the definition of the functions tfk and ipk (equations (|16[) and (|5ip) it is easy to 
rewrite the last four equalities as follows 



ScZ(">(Mo;z)= £ (-i)*(£W fc (a; 

even fc— 



)^fc(y) 
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and 

n 

h+1 n 



Sc Z^(j,z ;z)= J2 (-ir ! M^kM 



odd k=l 



Thus, for every n we have two nontrivial exact solutions of (|47[) (except for n — for which we observe 
that by construction Sc Z^(j, Zq;z) = 0). This infinite family of solutions can be written as follows 

u {x,y) = f{x)g(y), (58) 

m+l 

l 2 k / m+l \ 

u m {x,y) = ScZ [ ^ ) {l,z ;z) = ^2 ( £ ) fm^l- k ( x )^k(y) for an odd m, (59) 

oven fe=0 ^ ' 



2 fe+i /™\ 

Um (x,y)-ScZ^)(j,z ;z)= £ (-1)- ? 

-dd fe=l 



iprn.-k{x)ipk(y) for an even to. (60) 



Remark 24 In £/ie case when / = 1 and 5 = 1 we obtain that the system {u m }m=o * s ^ e system 
of harmonic polynomials {Sc(z — zo) n , Sc (J(z — ^o)")}^Lo - Theorems about its completeness like 
the Runge theorem and further related results are well known (see, e.g., J^f, J27jj, J30f . J31y). It is 
convenient to introduce the notation 

Po(x,y) = 1, 

m+l 

2 / m+l \ 

p m (x,y) = Sc{z - zo)^ = (-l)M j? J(^-^o) I22ti_fc (y-2/o) fc for an odd m, 

even k—0 

JL, fc+i /™\ „ 

Pm{x, y) = Sc (j(z - z ) 2 ) = 2^ (-1) 2 I 2 J (x - x ) 2 (y- y ) for an even m. 

odd fe=l ' ' 



Remark 25 Every u m is a result of application of an operator of transmutation to the corresponding 
harmonic polynomial p m . Indeed, consider for simplicity z = and suppose that f is defined on 
the segment [—a, a] and g is defined on [—6,6]. Both functions are assumed to be Ci-valued twice 
continuously differ entiable and nonvanishing. Let T f be the operator T defined by h28\) and T g be its 
equivalent associated with the function g. That is, 

T g v(y) = v(y) + f K g (y, t; g'(0))v(t)dt (61) 
J-v 

where 

K 9 (y, t; g'(0)) = ^ + K g ( y , t) + ^ £ {Kg ^ s) _ Kg(Vi _ s)) ds 
and Kg is a solution of the Goursat problem 

d 2 \ d 2 



1 



y 



K 9 (y,y) = 2 J o <ig( s ) ds > K g{y,-y) = 

with q g :— g" jg. Then 

u m (x, y) = T f T g p m (x, y). (62) 

This relation follows from Theorem \ll\ according to which the operator T/ maps a k-th power of x into 
ipk(x) for any k 6 No and similarly the operator T g maps a k-th power of y into ipk(y). Moreover, 

TfTgP m = T g Tfp m . 
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This observation together with the Runge approximation theorem for harmonic functions allows 
us to prove the following Runge-type theorem for the family of solutions {u TO }^ =0 . 

Theorem 26 Let O C R = [—a, a] x [—6, b] be a simply connected domain such that together with 
any point (x,y) belonging to Q the rectangle with the vertices (x,y), (—x,y), (x,—y) and (—x,—y) 
also belongs to £1. Let the equation 

(-A + q(x,y))u(x,y) = (63) 

in fl admit a particular solution of the form 4>(x,y) = f(x)g(y) where f and g are Ci-valued functions, 
f G C 2 [—a, a], g G C 2 [—6, b], f(x) ^ 0, g{y) ^ for any x G [—a, a] and y G [—b, b] (obviously, q has 
the form q(x,y) = qi(x) + 92(2/) with q\ = f" / f and q2 = g" /g). Then any solution u of UJ3j) in fl 
can be approximated arbitrarily closely on any compact subset K of Q by a finite linear combination 
of the functions u rn . That is for any e > there exists such a number M G N and such coefficients 

M 

< e for any point (x, y) G K . 



{ a m}m=0 C C i that 



u(x,y) - ^2 

(x,y) 



m=0 



-lm-l„ 
f L 9 

be approximated in K with respect to the maximum norm by a harmonic polynomial, \\v — Pm\\ < £1, 

M 

where Pm = ct m p m . Now we use the fact that T f and T g are bounded Volterra operators 

m— 

possessing bounded inverse operators. We have 

II" - T/T 9 P M || = \\T f T g v - T f T g P M \\ < £1 \\T f \\ \\T g \\ = e 

where the norms ||Ty|| and ||T ff || can be estimated in terms of the maximum values of the corre- 
sponding (continuous) kernels K/ and K s . ■ 

The restricting condition on the shape of the domain 17 is due to the necessity to have well defined 
the function v as the image of u under the application of the transmutation operators. Let one of 
the functions / or g be real-valued, for example, /. Then to prove the completeness of the system 
{ Um }m=o usm § a transmutation operator one can assume the symmetry of the domain only with 
respect to the variable y. 

Theorem 27 Let O C R 2 be a simply connected domain such that together with any point (x, y) 
belonging to the point (x,—y) and the segment joining {x,y) with (x,—y) belong to £1 as well. 
Let the equation i63\) in admit a particular solution of the form <p(x,y) = f(x)g(y) where f is a 
real-valued and g is a d-valued, both twice continuously differentiable and nonvanishing up to the 
boundary functions. Then any solution u of 163\) in £1 can be approximated arbitrarily closely on any 
compact subset K of £1 by a finite linear combination of the functions u m . 

Proof. Consider the equation 

(~A + qi (x))v(x,y) = inn (64) 

where qi = /"// is real valued. Denote by u m the functions defined by (|55)) - lpD|) with g = 1. Then 
as was proved in [5] any solution v of (|64p can be approximated arbitrarily closely on any compact 
subset of fl by linear combinations of the functions u m . As u m — T g u m , m G No (where g is the factor 
in 4> depending on y) once again using the boundedness of T ff and T" 1 we obtain the completeness 

Extension of the results of the preceding two theorems onto arbitrary simply connected domains 
is possible if equation (p3")) has the Runge property (see, e.g., [31], [1], [TU]). 
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Definition 28 Equation Lu = is said to have the Runge approximation property if, whenever, fii 
and 0,2 are two simply connected domains, f2i a subset of Q2, any solution in f2i can be approximated 
uniformly in compact subsets of fix by a sequence of solutions which can be extended as solutions to 

n 2 . 

It is known [4], [10] that the Runge property in the case of elliptic equations with real- valued 
coefficients is equivalent to the (weak) unique continuation property (if every solution of Lu = which 
vanishes in an open set vanishes identically) and is true, e.g., for second-order elliptic equations with 
real-analytic coefficients. Without going into further details concerning the Runge approximation 
property which is beyond the scope of the present work, we prove that if equation (|63p has this 
property then the family of solutions {a m }™ =0 is complete in any simply connected domain. 

Theorem 29 Let equation A63\) in a rectangle R = (—a, a) x (—6, b) admit a particular solution of the 
form <j)(x,y) = f(x)g(y) where f and g are arbitrary Ci-valued twice continuously differentiable and 
nonvanishing functions in [—a, a] and [—b, b] respectively. Let £1 C R be a simply connected domain. 
Assume equation 163\) has the Runge property. Then any solution u of i63\) in Q can be approximated 
arbitrarily closely on any compact subset K of Q by a finite linear combination of the functions u m 
defined by if55|)-fg^. 

Proof. Consider a solution u in Q which due to the Runge property can be approximated on if by a 
solution v of (|63[) in R. Due to Theorem 1261 v in its turn can be approximated on K by the functions 
u m from where we obtain the required approximation of the solution u in terms of the solutions u m . 



In the rest of the present section we show that solutions of (|63j) sufficiently smooth up to the 
boundary of the domain of interest Q can be approximated by functions u rn in f2. By we denote 
the linear space of solutions of (|63|) in il satisfying the following regularity requirement u € C 2 (f2) n 
C 1+Q (£l), < a < 1. This linear space can be equipped with one of the following scalar products 
(we assume that zero is neither a Dirichlet nor a Neumann eigenvalue) 

f f du dv* 

<u,v>i= / uv*ds and < u, v >2= It, — ^ — ds (65) 

on on 



an 



where by "*" we denote the complex conjugation in Q, and -g- is the outer normal derivative. With 
the aid of the scalar products (|65[) two Bergman-type reproducing kernels [2] can be introduced for 
solving the Dirichlet and Neumann problems respectively as well as the corresponding eigenvalue 
problems [5]. 

A complete orthonormal system of functions in with respect to < •, • >i or < •, • >2 allows 
one to construct a corresponding "Dirichlet" or "Neumann" reproducing kernel respectively. In 
[5] the completeness of the family of solutions {u m }^ =0 of (|63[) obtained as real parts of complex 
pseudoanalytic formal powers was proved in the case when (163[) admits a particular solution in a 
separable form (f{s,£j — S(s)T(t) where S and T are arbitrary twice continuously differentiable 
nonvanishing real- valued functions, $ = s+it is a conformal mapping defined in and O is a domain 
bounded by a Jordan curve. In the same paper it was shown that the completeness of {u m }^ =0 in 
Tfl a in the case when the particular solution ip is complex- valued is an important open problem and 
its solution is required not only for solving boundary value problems for (163[) with a complex- valued 
coefficient but also for solving spectral problems for (|63]) even in the situation when the coefficient 
is real-valued. Here by means of the developed results concerning the transmutation operators we 
obtain the completeness of the family of solutions {« m }™ =(J in E^ under the conditions of Theorems 
HUandHZl 

Theorem 30 Let O C R = [—a, a] x [—6,6] be a simply connected domain such that together with 
any point (x,y) belonging to f2 the rectangle with the vertices (x,y), (—x,y), (x,—y) and (— x, — y) 
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also belongs to O. Let equation \6S\) in admit a particular solution of the form <f>(x,y) = f{x)g{y) 
where f and g are Ci-valued functions, f € C 2 [—a, a], g G C 2 [—b,b], f{x) ^ 0, g(y) ^ for any 
x G [—a, a] and y G [—b,b]. Then the family of solutions {w m }^ =0 is complete in E^ 7 a > tuif/i 
respect to both norms generated by the scalar products \65\) . 

Proof. Let u G E*. Consider the harmonic function v := 1 T~ 1 u which belongs to E° due to the 
fact that the kernels in both transmutation operators are at least C^-functions and hence the operator 
T^T" 1 transforms E' into E°. There exists (see, e.g., 0) a sequence of harmonic polynomials Pm 
such that when M — > oo, Pm — > v uniformly in Q, together with their first partial derivatives. It is easy 
to see that this implies the uniform convergence in VL of the sequences T/T s Pa/ = Um — > u = T/T g w, 
—> and — > Indeed, the uniform convergence of Um to u follows directly from the 
boundedness of the transmutation operators, see the proof of Theorem 1261 and the verification of the 
uniform convergence of the partial derivatives is straightforward. Consider 

—T f P M {x,y) = —P M (x,y) + —J K f (x,t;h)P M (t,y)dt 
d [ x d 

= T^PM(x,y) + J T^K f (x,t- 1 h)PM(t,y)dt + Kf(x,x;h)P M {x,y) + K f (x, -x; h)P M (-x, y) 

from where due to the uniform convergence of Pm to v together with their partial derivatives and 
due to the fact that -j^K.f{x, t; h) is continuous, it follows that ^j^f- — > uniformly in f2. The proof 
of the uniform convergence of S^m. \ ^± j s analogous. 

Now, the completeness of {« m }™ =0 in E* with respect to the norm generated by < •, • >i follows 
from the uniform convergence of Um to u in Q, and hence from the completeness of {w m }^ =0 with 
respect to the maximum norm in fl. To verify the completeness of {w m }^ =0 in E^, with respect to 
the norm generated by < •, • >2 consider the following chain of relations 

\\u-U M \\= f d ( u - U M)d(u*~U* M ) ds ^ r | V(u _r7 M) |2 ds < isup |v( u -tf M )| 2 ^ . 
Jan On on J gn dn 



Theorem 31 Let O C R 2 be a simply connected domain such that together with any point (x, y) 
belonging to D, the point (x,—y) and the segment joining {x,y) with (x,—y) belong to as well. Let 
equation i63\) in Q admit a particular solution of the form <f>{x,y) — f{x)g(y) where f is a real-valued 
and g is a Ci-valued, both twice continuously differentiable and nonvanishing up to the boundary 
functions. Then {u m }^ =0 is complete in E*, a > with respect to both norms generated by the 
scalar products A65\) . 



Proof. The first part of the proof is similar to that of Theorem \T7\ We have that u m = T g u m , 
m G No where {u m }™ =0 is complete in EJ? 1 (see [5]) with respect to the required norms. Then the 
completeness of {u m }™ =0 is proved analogously to the proof of Theorem [301 ■ 

Remark 32 In order to prove the completeness of the family of solutions {w m }^ =0 in E-j, with respect 
to both norms generated by the scalar products \65\) under less restrictive conditions on the shape of 
the domain Q in fact we need a result on the existence of an infinite system of solutions of i63\) in 
R D 51 and complete in E* (<9f2) or in a maximum norm inQ. If such a system exists then according 
to Theorem \26\ every element of it can be approximated arbitrarily closely by linear combinations of 
functions u m which would allow one to prove the completeness o/{it m }^ =0 . Thus, if such a complete 
system exists then {« m }™ =0 is precisely such system. The question on the existence requires further 
study. 
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7 Conclusions 



Transmutation operators for Sturm-Liouville equations are considered and their new properties con- 
cerning the transformation of certain infinite systems of functions generated by the Sturm-Liouville 
operators are presented. These infinite systems of functions slightly generalize the notion of L- 
bases [TT] and play an important role in the theory of linear differential equations. We show how a 
transmutation operator can be constructed mapping one such basis into another and give an appli- 
cation of this result obtaining several theorems on the completeness of certain families of solutions 
of two-dimensional stationary Schrodinger equations which are obtained as scalar parts of bicomplex 
pscudoanalytic formal powers. To our best knowledge this is the first result of this kind in bicomplex 
pseudoanalytic function theory. Its importance is in the fact that it opens the way for construction 
of Bergman-type reproducing kernels for corresponding second-order elliptic equations with variable 
complex- valued coefficients and hence for solving boundary and eigenvalue problems. 
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